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We have completed a new measurement of the Stark shift in ^^^In within the 410 nm 5Pi/2 — > 
6S1/2 transition. We measure the Stark shift constant to be ks ~ —122.92(33) kHz/(kV/cm)^, 
corresponding to a difference in the 6S1/2 and 5Pi/2 state polarizabilities, Aqq, of 1000.2 ± 2.7 Uq 
(in atomic units). This result is a factor of 30 more precise than previous measurements and 
is in excellent agreement with a new theoretical value based on an ab initio calculation of the 
wave functions in this three- valence-electron system. The measurement was performed in an indium 
atomic beam apparatus, used a GaN laser diode system, and exploited an FM spectroscopy technique 
to extract laser transmission spectra under conditions where our interaction region optical depth 
was typically less than 10~^. 



I. INTRODUCTION 

High precision atomic structure measurements have 
long been an important tool in testing the accuracy and 
guiding the refinement of theoretical techniques aimed 
at calculating atomic wave functions in multi-electron 
atomic systems. Alkali systems have played a partic- 
ularly important role in this experiment-theory inter- 
play. More recently, new ab initio calculational tech- 
niques have produced theoretical results of improved pre- 
cision for multi- valence-electron systems [IHl]. Of rele- 
vance here arc the Group III systems such as indium 
and thallium. Thallium, in particular, has played an im- 
portant role in tests of discrete symmetry violation over 
recent decades[5HZ], and a proposal to measure the per- 
manent electric-dipole moment of indium was recently 
published JS]. The size of these symmetry- violating ob- 
servables scales rapidly with the atomic number, encour- 
aging the use of high-Z systems. This therefore requires 
independent, precise atomic wave function calculations in 
order to distinguish quantum mechanical effects from the 
elementary particle physics observables being targeted. 
For example, the present 2-3% uncertainties in ab ini- 
tio wave function calculations in thallium currently limit 
the quality of the standard model test provided by a 1995 
thallium parity nonconservation measurement fB' . A very 
similar theoretical method to that used for thallium can 
be applied to other three- valence systems such as indium 
and gallium [2]. 

Recently our group measured the hyperfine constants 
of the indium 6P3/2 excited statej9 , providing a test of 
short-range electron wave function behavior as well as a 
measurement of the nuclear quadrupole moment. The 
Stark shift result reported here is quite complementary 
to this hyperfine structure work in that it provides a test 
of long-range wave function properties. Previous exper- 
imental work to measure indium polarizability includes 
a 1970 atomic beam spectroscopy experiment [10]. and a 



* pmajumdeOwilliams.edul 



1984 ground state measurement based on atomic beam 
deflect ion [TT]. both at the 10% level of accuracy. 

Our new polarizability measurement, with its 0.3% ac- 
curacy provides a benchmark test of two distinct, high- 
precision calculational methods which can be applied to 
multi-valence atoms such as indium. One method, em- 
ploying a coupled-cluster approach, has been particularly 
effective for mono-valent systems [T2], but also produced 
results l3j in excellent agreement with experimental polar- 
izability measurements in thallium |13l[T^ . A second the- 
oretical strategy incorporates a configuration interaction 
approach to better accommodate three-particle states. 
Both of these methods were very recently employed to 
compute the 6S1/2 and 5Pi/2 state polarizabilities in in- 
dium. The agreement between these theoretical methods 
and between theory and our experimental result is very 
good, and is detailed in [15]. Finally, as discussed at the 
end of this manuscript, our experimental result and the 
theoretical expressions can be combined to extract new, 
accurate values for the indium 6P1/2 and 6P3/2-state life- 
times. 



II. ATOMIC STRUCTURE DETAILS 

We focus on the principal naturally-occuring isotope 
of indium (^^^In, 96% abundant), which in our Doppler- 
narrowed atomic beam geometry is spectroscopically iso- 
lated from the small ^^^In component. ^^^In has nu- 
clear spin / = 9/2, so that both the 5Pi/2 and 65*1/2 
states contain F = 5 and F — 4 hyperfine levels. The 
respective 11.4 and 8.4 GHz hyperfine splittings (HFS) 
of the ground and excited states again yield an entirely 
resolved spectrum in our atomic beam apparatus. Be- 
cause we study aJ = 1/2— J = 1/2 transition, the 
Stark shift of each level has only a scalar component, 
producing a common shift of all sublevels within a given 
state, and yielding an experimental result which is inde- 
pendent of relative laser and static electric field polar- 
ization. Expressing the energy shift of a given level as 
AW = —^oqE'^, where ag is the scalar polarizability. 
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the observed frequency shift of the 410 nm hne can then 
be expressed as: Ai/g = — 57^ [ao (651/2) — cto{5Pi/2)]E^ ■ 
In the one electron central field (OECF) valence approxi- 
mation, the polarizability is calculable from second-order 
perturbation theory, and for the state 'v' is given by the 
expression [T6| [T7] 
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where (fc||Z?||u) is the reduced electric-dipole (El) matrix 
element, and the sum extends over all states for which 
the El matrix elements are non-zero. 

Due to the size of the HFS compared to the optical 
transition frequencies {Svhfs /^^opt of order 10^^) we ex- 
pect no measurable difference in Stark shift value among 
four hyperfine component lines, given our level of exper- 
imental precision. Small tensor components to the Stark 
shift of individual sublevels induced by higher-order hy- 
perfine interaction effects are similarly of negligible mag- 
nitude in our experiment fTS^. As discussed below, we did 
acquire data for two different hyperfine transitions, and 
different laser linear polarization directions, as part of 
our overall search for potential systematic errors. 

In this experiment, we perform transverse laser spec- 
troscopy. Given the geometry of our atomic beam and in- 
teraction region, we find that the residual Doppler broad- 
ening of our atomic spectrum is reduced roughly 15-fold 
from the value that would be found in our heated source. 
As such, the absorption line shape is described by a Voigt 
profile with a 23 MHz Lorentzian component due to the 
natural line width, and a roughly 100 MHz Gaussian 
component due to the residual Doppler broadening. As 
described below, we use an electro-optic modulator to 
add FM sidebands to the spectrum of our blue laser. We 
demodulate the subsequent atomic beam transmission 
signal at either the first or second harmonic of the modu- 
lation frequency. The resultant atomic line shapes can be 
computed as a function of FM modulation index, choice 
of demodulation frequency, and relative phase [19l [20] . 
As discussed in Sec. IV, we have developed a model for 
exact fitting of our RF-demodulated spectra, but also 
perform simplified polynomial or Lorentzian fits to the 
central portion of our line shapes to extract line centers 
and Stark shift information. 



III. EXPERIMENTAL DETAILS 
A. Atomic beam system 

Fig. 1 shows a cross section of our atomic beam oven 
source, vacuum chamber, and interaction region. Our 
atomic beam source consists of a cylindrical molybde- 
num crucible capable of containing 100 grams of metal- 
lic indium. The crucible is supported by two legs that 
tilt the cylinder at roughly 20° upwards from the hori- 
zontal. The open front face of the cylinder is milled so 




FIG. 1. A cross-sectional view of the atomic beam source, vac- 
uum chamber, and interaction region. The source-interaction 
distance is roughly 40 cm. 



that it is in the vertical plane, and we cap the crucible 
with a molybdenum faceplate containing 20 parallel 0.25- 
mm-wide vertical slits extending over a 2 cm horizontal 
width. The 0.5 cm thickness of the faceplate serves to 
help pre-coUimate our atomic beam by creating a series 
of parallel 'tunnels' from which the indium atoms must 
effuse. The crucible sits within an alumina tube around 
which we attach a pair of resistive heaters. These are 
surrounded by layers of carbon felt insulation and thin 
stainless steel (serving as a heat shield). A water-cooled 
cylindrical copper cold wall forms the outermost layer of 
our source oven. The whole assembly rests on an ad- 
justable base that bolts to the floor of the vacuum cham- 
ber. Four thermocouples attached to crucible, heater, 
stainless steel heat-shield and copper cold-wall allow for 
continuous temperature monitoring. To reach our tar- 
get temperatures between 900 and 1000°C, we require 
roughly 0.6 kW of (AC) power. We regulate the tem- 
perature of the crucible using a variable-duty-cycle con- 
troller which can produce stable source temperatures up 
to 1100°C. 

Our stainless steel vacuum chamber includes a pair of 
liquid-N2-trapped diffusion pumps (Varian M-6), one lo- 
cated beneath the source oven, and one beneath the in- 
teraction region chamber. Under hot-source-oven oper- 
ating conditions, we typically observe lxlO~^ torr pres- 
sure in our source chamber, while the pressure in the 
differentially-pumped interaction region is 3 or 4 times 
lower. Roughly half way along the ~ 40 cm source- 
interaction region path we place adjustable razor blade 
collimators to define and coUimate our atomic beam in 
both vertical and horizontal directions, creating a 2-mm- 
high, 1.5-cm-wide ribbon-shaped atomic beam. Just up- 
stream of our electric field plates we insert one additional 
coUimation iris of comparable dimensions to minimize the 
amount of scattered metallic indium which might coat 
the field plate surfaces. Finally, we attach an in-vacuum 
chopping wheel to the final collimator assembly, just up- 
stream of our interaction region. In this way, we can mod- 
ulate the atomic beam at a typical frequency of 500 Hz, 
leading to removal of background features in our atomic 
spectra as discussed below. A set of three pairs of or- 
thogonal magnetic-field-cancellation coils on the outside 
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of the interaction chamber is used to reduce the ambient 
magnetic field to < 10~^ T within the interaction region. 



B. Electric field plates and high voltage system 

Our electric field plate assembly consists of a pair of 
highly polished circular stainless steel plates that are 10 
cm in diameter. The two plates are separated by four 
ceramic spacers and positioned horizontally on a ceramic 
pedestal. The separation between the plates was repeat- 
edly measured using metric gauge blocks to be 0.999(1) 
cm. The modest temperature increase characteristic of 
the downstream interaction region during beam opera- 
tion causes negligible change in plate separation, espe- 
cially given that expansion of the ceramic posts and steel 
plates tend to oppose one another. Our blue laser beam 
intersects the 1.5-cm-widc atomic beam at the center of 
the plate assembly. Independently, we have performed 
simulations to insure that, given our interaction region 
geometry, any non-idealities due to the finite size of the 
plates would affect our electric field calibration only at 
the level of 1 part in 10^ or below. 

A high voltage power supply (Glassman ER40P07.5) 
allows the application of up to 40 kV to the plates. Our 
data acquisition system allows computer control of both 
polarity and magnitude of this voltage. In series with the 
HV output and the field plates, we place both a chain 
of ballast resistors of total resistance 100 MJ7 (for cur- 
rent limiting and power supply protection), as well as a 
1 MCI resistor between the ground plate and the power 
supply retTirn (for leakage current monitoring). Leakage 
currents were typically in the nA range, so that errors 
due to leakage current-induced voltage drops were negli- 
gible for our purposes. In parallel with the power supply 
we installed a commercial high voltage divider (Ross En- 
gineering, Inc.). This device has a quoted precision of 
0.01%, a total resistance of 180 MQ, and is customized 
to be impedance matched to the particular voltmeter we 
use (Keithlcy model 197A, 0.011% absolute voltage mea- 
surement accuracy). The data acquisition program reads 
the output of the Keithley voltmeter via a GPIB inter- 
face, and, to the quoted precision of the instruments, this 
is exactly 1/1000 of the actual high voltage applied. 

C. Optical System 

The experiment uses a commercial external cavity 
GaN diode laser system in the Littrow configuration 
(Toptica Photonics, DL 100) to produce roughly 10 mW 
of light at a wavelength near 410.2 nm. The laser can be 
tuned to excite any of the four hyperfine transitions in 
indium 5Pi/2(F = 4, 5) ^ 6S1/2 (F' = 4, 5) manifold. 
After passing through an optical isolator, a fraction of 
the laser output is directed into a Doppler-free saturated 
absorption system. This system includes an optical 
chopping wheel which interrupts the pump beam, and 



we record the subsequent lock-in amplifier output of the 
probe beam signal. We do not lock our laser using this 
signal, but instead use the Doppler-free absorption peak 
as a frequency reference to facilitate comparison of the 
spectral shift in atomic beam transmission signal when 
we turn on and off the electric field in our interaction 
region. 

The main laser beam is then passed through an electro- 
optic modulator (EOM; New Focus, Model 4001). The 
EOM is driven by a 100 MHz RF synthesizer and power 
amplifier. We obtain a frequency modulated (FM) laser 
spectrum whose exact sideband distribution is deter- 
mined by the EOM input power. A portion of the 
modulated laser beam is then directed to a 1 GHz free- 
spectral-range confocal Fabry-Pcrot (FP) Interferometer 
(Burleigh RC-110) whoso finesse is roughly 40. We isolate 
and passively temperature stabilize the FP cavity by cov- 
ering it with acoustical insulation material. By detect- 
ing the FP transmission signal as wo scan our laser fre- 
quency, the FM laser spectrum, with its precisely-known 
100 MHz sidebands, can thus be monitored, and the laser 
frequency scan can be calibrated and linearized. In order 
to create a spectrum with measurable second and third- 
order FM sidebands, we drive the EOM with sufficient 
power to achieve modulation depths, /3, ranging from ap- 
proximately 1.0 to 3.0. 

The modulated laser beam is then directed into and 
through our atomic beam apparatus via windows which 
are anti-reflection coated for our 410 nm wavelength. As 
it enters the vacuum chamber, the laser beam has a 1- 
mm diameter, an adjustable linear polarization, and is ar- 
ranged to intersect the atomic beam in a transverse fash- 
ion. By passing our laser beam through 1-mm-diamcter 
collimating apertures on either side of our vacuum cham- 
ber, we insure that the we have a reproducible, well- 
defined interaction geometry. This minimizes any pos- 
sible Doppler-shift-related systematic errors in our spec- 
tra. After exiting the vacuum chamber, having interacted 
with the atomic beam, we collect the optical transmission 
signal on a 1-GHz-bandwidth photodiode (New Focus, 
model 1601). Fig. 2 shows a complete schematic layout 
of our optical and signal processing arrangement. 



D. FM signal processing and lock-in detection 

By applying a linear voltage ramp to the piezoelectric 
transducer (PZT) of the diode laser external cavity, we 
scan the 410 nm laser frequency upwards and downwards 
over ^1.5 GHz which includes either the F = 5 F' = 
4orF = 4— >F' = 5 hyperfine resonance. Because the 
vapor pressure of indium is comparatively low (for exam- 
ple, 1000 times lower than thallium in this temperature 
range) , we observe an atomic beam optical depth of ap- 
proximately 0.001 in our interaction region even at oven 
temperatures near 1000°C. We have therefore made use 
of an FM spectroscopy / RF detection technique which 
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FIG. 2. A complete experimental setup indicating the config- 
uration of the Doppler-free saturated absorption system (top), 
Fabry-Perot interferometer, and transverse atomic beam FM 
spectroscopy system (below). 



greatly improves the quality of our atomic absorption 
spectra under conditions where direct transmission mea- 
surement is not a feasible option. As described below, we 
demodulate the transmitted atomic absorption signal at 
either the fundamental FM frequency ('If = 100 MHz 
in our case) or at the second-harmonic frequency (2f). 
One can also study the in-phase or the quadrature com- 
ponent of the demodulated signal. We find that higher 
FM modulation depth complicates the line shape, but 
provides us with a larger signal size. A complete analysis 
of the form of the expected demodulated signals under 
various conditions can be found elsewhere [THl I20j. The 
frequency-modulated electric field can in general be ex- 
pressed in terms of sideband components with relative 
amplitudes given by Bessel functions, J„, whose argu- 
ments are the modulation depth. Defining and Wm 
to be the carrier and modulation frequency respectively, 
and (3 to be the modulation depth, we have: 
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After interacting with the atomic beam, we can as- 
sociate a (complex) electric field transmission function, 
Tn{i^c + nuJrn), with cach sideband component. In the 
limit of very small absorption, we approximate this as: 
Tn ~ I — {Sn + i4>n), whcrc S is the frequency-dependent 
atomic absorption line shape and is the associated (dis- 
persive) phase shift. The transmission signal we detect is 



then proportional to the absolute square of the transmit- 
ted field. Our RF lock-in amplifier (SRS model SR844) 
is referenced to the synthesizer which drives the EOM. 
By detecting the Fourier component of the transmitted 
intensity at either If or 2f, we pick out particular cross 
terms in the square of the sum above. We can derive an- 
alytic expressions for the in-phase and quadrature com- 
ponent demodulated spectra of both the If and 2f de- 
modulated signals for arbitrary modulation depth. For 
example, for the case of demodulation at the fundamen- 
tal modulation frequency, we find [20]: 

oo 
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where 9d is the relative phase between the FM signal 
and the lock-in detector. Each absorption component 
(5„(a;c + rujJm) represents a distinct, frequency-shifted 
Voigt profile. We use an appropriately truncated sum of 
this form to fit our demodulated spectra, where the mod- 
ulation depth and lock-in phase are left as fit parameters 
in our analysis. 

Eq. ^ predicts a zero-background demodulated sig- 
nal. Because the absorption in our atomic beam is so 
small, in practice we find that our demodulated spectra 
contain some non-zero frequency-dependent background 
in addition to our atomic absorption signal, likely due 
to some purely optical effects, such as etaloning due to 
small reflections within the EOM itself. We can com- 
pletely eliminate this residual optical background pat- 
tern using our atomic beam chopping wheel. We direct 
the output signal of our RF lock-in amplifier to a sec- 
ond low-frequency lock-in amplifier (SRS model SR810) 
which is referenced to the chopping wheel rotation fre- 
quency. The output of this second-stage demodulation is 
a truly background-free line shape which we record and 
analyze. 



E. Data acquisition and experiment control 

The high voltage (HV) is programmed and controlled 
using a Lab VIEW program. A computer-generated tim- 
ing sequence initiates switching on or off the high voltage. 
After each state switch, we wait roughly 14 seconds for 
all high-voltage-switching transients to subside. This de- 
lay time was chosen by studying the consistency of scans 
acquired with various shorter delay times. We then ini- 
tiate a triangular external voltage ramp signal which we 
apply to the PZT in the laser's external cavity. Two up- 
ward and downward going sweeps, each requiring roughly 
3 seconds, are completed for each high voltage state. In 
this way, we can compute the Stark shift by comparing 
consecutive laser scans in two different ways. By compar- 
ing the [HV on - HV off] frequency shift to that obtained 
with the temporal order reversed we can study and po- 
tentially eliminate any slow drift, or HV state-specific 



5 



C D 



Laser Ramp 
Data Acq. 



FIG. 3. A timing schematic for the experiment and data 
acquisition. We compute a set of Stark shift values by com- 
paring HV on-off scans, such as B-C above, and a separate 
set of values by comparing HV off-on scans such as D-E in 
the figure. 



effect in the experimental system that might lead to a 
systematic error in our measured frequency shift. In all 
cases, we record and analyze upward- and downward- 
going laser scans separately. Fig. |3] summarizes the data 
acquisition timing sequence. The FM signal, FP trans- 
mission signal, Doppler free saturated absorption (SA) 
signal, and the output of the meter which monitors the 
field plate high voltage are sent to our data acquisition 
interface and recorded by the Lab VIEW program. 
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FIG. 4. Simultaneous data acquisition of demodulated FM 
signal (top), Fabry- Perot transmission, and saturated absorp- 
tion signal (bottom). For this scan, we demodulate the FM 
signal at If, and the entire scan extends over roughly 1.5 GHz. 



B. Data fitting procedure 



IV. DATA AND ANALYSIS 
A. Summary of data collected 

We typically acquired alternating HV off / HV on 
scans in the sequence described above using a particular 
choice of HV for a 30-minute period. The actual volt- 
age measured by our high-precision Keithley voltmeter 
is recorded for each data point, though in practice the 
voltage value varies only in the third decimal place from 
cycle to cycle within a set of nominally identical HV 
runs. This produced ^ 200 individual spectra, provid- 
ing ^ 100 Stark shift determinations, of which half are 
upward-going laser scans, and half downward scans. We 
then changed our nominal value for the HV, repeating 
these half hour data collection cycles for 11 values of HV 
between 10 kV and 20 kV. We analyzed 18 such days of 
data in total. A wide variety of experimental conditions 
were changed over the course of the full data-collection 
period, including hyperfine transition, laser polarization 
and power, EOM power (modulation depth for FM spec- 
trum), and RF lock-in detection frequency (If, 2f). In 
total, we obtained roughly 6,000 individual Stark shift 
measurements which were analyzed in multiple ways to 
insure that the Stark shift values extracted from the fits 
were consistent and reliable. The final statistical error in 
our measurement is less than 0.1%. The high statistical 
precision thus allowed us to explore many possible sys- 
tematic error contributions to the data set, as we discuss 
below. 



Since in this experiment we compare atomic resonance 
positions after switching HV state, we require that the 
system remain stable over time scales of at least tens of 
seconds. Furthermore, we require a spectral feature that 
can be used as a frequency reference point. Finally, as 
we ramp the voltage applied to the laser cavity PZT for 
our scans, we must linearize and calibrate our scan to 
produce an accurate frequency axis. For each laser scan, 
we simultaneously collect the demodulated FM atomic 
beam signal, the FP transmission signal, and the sat- 
urated absorption signal from our supplementary vapor 
cell. Fig. |4] shows such a single scan and the associated 
signals that we collect. For this data scan we employed If 
demodulation with an FM modulation depth of roughly 
1.5. We typically scan the laser over a ^ 1.5 GHz range 
centered on the atomic resonance. To begin a data col- 
lection cycle, we manually tune the FP cavity length to 
insure that the FM spectrum of the FP cavity is centered 
on the atomic resonance, and we find that we do not need 
to retune the cavity for several hours. We typically col- 
lect 400 data points per 3-4 second scan. Prior to further 
analysis, we recast our frequency axis in terms of a nor- 
malized point number, xj = [j — N/2) /{N/2). In this 
case N = 400, and 1 <j < 400, so that -1 < Xj < +1. 

We begin the analysis procedure by fitting the cen- 
tral FP spectrum. Because we limit our laser scan to a 
single set of FM sidebands associated with a single FP 
longitudinal mode, we need not fit our FP spectrum to 
a periodic Airy function and instead fit the spectrum to 
a simple sum of seven Lorenztians. We know that each 
Lorenztian peak is separated from the next by precisely 
100 MHz, our modulation frequency. However, due to 
nonlinearity in the PZT response, we expect that the fit- 
ted Lorenztian peak locations, x^^-* . . .x^^-*, will not have 
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equal separation in terms of the variable x. Given the rel- 
atively small frequency-scanning range, and the nonlin- 
earity typical for such PZT devices, we find that a third- 
order polynomial allows an accurate frequency mapping 
of our scan. Thus, to each point in our scan we assign a 
frequency according to /(j) = oq + aiXj + 022;^ +030;^, 
where the polynomial coefficients are determined by in- 
sisting that the Lorenztian FP peaks occur at exact 100 
MHz intervals. Fig. [5] shows the results of a fit to an 
FP spectrum. The deviation from linearity in our scan is 
small, but we find that quadratic and cubic components 
can be resolved. These parameters tend to remain consis- 
tent for repeated scans taken under identical conditions. 
Moreover, upward- and downward- going laser scans tend 
to have opposite signed quadratic components due to the 
hysteretic nature of the PZT. We found that for these 
scans, a fourth order polynomial fit did not reliably re- 
solve a quartic coefficient and did not improve the quality 
of the fit. While this polynomial model may not be a reli- 
able mapping of the frequency near the edges of the scan, 
where FP peaks do not exist, it does provide an accurate 
map of the frequency non-linearity in the central portion 
of the scan, exactly where the FM atomic beam signal is 
located. 
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FIG. 5. A Fabry-Perot transmission data (points) with the 
results of a fit to a sum of seven Lorentzians (solid line). The 
seven peak locations are exactly 100 MHz apart. All of these 
peaks correspond to a single FP longitudinal mode (free spec- 
tral range = 1.0 GHz). 

In the ~14 seconds between HV on and off scans, a 
small change in the laser cavity PZT position or FP cav- 
ity length would obviously result in a systematic error in 
a Stark shift determination. For this reason, we use the 
saturated absorption (SA) peak from our supplementary 
vapor cell as a frequency reference for both scans. Using 
the newly calibrated frequency axis, we fit the central 
portion of the SA peak to a Lorenztian to determine the 
line center. Independently, we can evaluate the extent of 
possible slow drifts in the FP or PZT systems by com- 
paring the SA line center position to Oq (derived from 
the FP fit) over time scales up to one hour. We find that 



the maximum drift in relative line positions of the two 
signals to be 10 MHz over these time scales, and typically 
much less than this. Nevertheless, as the Stark shifts we 
measure are of order 10-50 MHz, the use of the SA cell 
reference is clearly essential. After the SA fit, we reset 
the frequency offset ao so that the SA peak frequency de- 
fines the frequency origin, while keeping the higher order 
parameters ai — 03 intact. We can now compare the line 
center of our FM atomic beam spectra to the location of 
the SA peak using a properly linearized and calibrated 
frequency scale. 

The Stark shift is defined as the difference between the 
fitted center of the demodulated atomic beam transmis- 
sion spectrum with the high voltage on and with the high 
voltage off (we observe a negative frequency shift upon 
application of the field) . If we demodulate the FM sig- 
nal at If (2f) the resultant line shape is an odd (even) 
function. At low modulation depths, the If signal has 
a simple dispersion shape, while the 2f signal consists of 
a symmetric peak with small additional lobes on each 
side of resonance. For larger modulation depths, the de- 
modulated signals grows in size, but also becomes sub- 
stantially more complicated. We have developed a curve- 
fitting algorithm based on functions derived from the full 
FM analysis outlined above. That is, using Eq. (§ as 
a basis for the fitting model, and truncating the infinite 
sum appropriately, we fit our demodulated data to a sum 
of symmetric and dispersive Voigt profiles, allowing the 
modulation depth, the relative phase of lock-in detection, 
and the Voigt profile component widths to become fit pa- 
rameters in a non-linear least squares fitting routine. Ex- 
amples of these 'full FM model' fits to our demodulated 
line shapes are shown in Figs. |6] and [7] For both the 
If and 2f-demodulated cases, we show a pair of spectra 
with and without a 20 kV/cm electric field applied. 

While these full FM fits do an excellent job of captur- 
ing the complexities of the line shape, we also pursued a 
simpler, computationally faster fitting algorithm in order 
to assess the robustness of the line shape fitting results. 
At low modulation depths, we can fit the central por- 
tions of the If and 2f line shapes to a linear function or 
a Lorenztian function respectively. For somewhat larger 
values of /3 one can develop a correspondingly more in- 
volved fitting function by fitting the central portion of the 
line shape to a higher-order polynomial, an odd function 
for the If signals and an even function for the 2f signals. 
Since, as we observe, the pair of Stark-shifted line shapes 
are identical in size and shape, differing only in frequency 
offset, it is reasonable to assume that any imperfection in 
the polynomial approximation approach will affect each 
spectrum identically and thus not lead to any system- 
atic error in the Stark shift determination. Furthermore, 
we can compare our results for If and 2f demodulation, 
and for low and high modulation depth to insure that 
Stark shift values we extract are not sensitive to the ex- 
act shape of the resonance line shape. By exploring the 
degree of the polynomial fit, and the extent over which we 
fit the central portion of the line shape, we insured that 
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FIG. 6. Two up scan FM signals demodulated at If. Trian- 
gles(circles) represent data taken with HV off (on). Solid lines 
are fits based on the theoretical model discussed in the text. 
Expanded residuals are shown below. 
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FIG. 7. Two up scan FM signals demodulated at 2f. Trian- 
gles(circles) represent data taken with HV off (on). Solid lines 
are fits based on the theoretical model discussed in the text. 
Expanded residuals are shown below. 



the final Stark shift results were independent of these de- 
tails. The line shape centers and Stark shift constants 
extracted from the two different fitting methods were in 
excellent agreement. 



C. Analysis of Stark shift results 

A typical single Stark shift measurement, determined 
from a pair of consecutive scans, yields a 3-5 MHz uncer- 
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Electric Field Squared (kV/cmf 

FIG. 8. Measured Stark shift vs. the square of the electric 
field. The result of one day's worth of data is shown, with 
all runs taken at each of 10 different nominal electric field 
values averaged. Statistical errors of individual points are 
comparable to the point size. The Stark shift constant, ks, 
(slope of fitted line) has an uncertainty of roughly 0.4% for 
these data. 



tainty, while the magnitude of the shift varies between 
10 and 50 MHz for the HV values we use. We store re- 
sults for upward laser sweeps ('up scans') and downward 
sweeps ('down scans') separately. We further separate 
our Stark shift results into data sets based on the order 
of field switching - that is, we have separate sets of on-off 
as well as off-on Stark shift results for each HV setting. 
For each experimental condition, we compute a weighted 
average and standard error for the measured Stark shifts 
over one 30 minute cycle at a single HV value. By di- 
viding the resulting Stark shift value by the square of 
the electric field for this data set we can obtain a re- 
sult for the Stark shift constant ks = Av/E'^, in units 
of kHz/(kV/cm)^, which can eventually be converted to 
the atomic scalar polarizability. Alternatively, after an- 
alyzing similar data sets for a full range of high voltages 
between 10 and 20 kV, we can plot our measured shifts 
vs. the square of the electric field, as is done in Fig. [8] for 
one complete day's worth of data. It is encouraging to 
see the expected linear relationship, and that the linear 
fit yields an intercept value statistically consistent with 
zero. When we extract a value for the Stark shift con- 
stant independently for each individual 30 minute cycle, 
for each value of high voltage, and under all experimental 
conditions, we obtain the histogram shown in Fig. |9]con- 
sisting of roughly 500 distinct ks values. We can obtain a 
final statistical mean and error from our data in several 
ways - for example, by analyzing individual data such 
as in Fig. [9j or by computing slopes such as indicated 
in Fig. |8] for each day, and then averaging over all 18 
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days of data collection. These various statistical analysis 
schemes gave consistent final results. This then leads us 
to a final mean and one-standard-devation statistical er- 
ror for our Stark shift measurement of: ks = —122.92(05) 
kHz/(kV/cm)2. 
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FIG. 9. A histogram of all ks values. Each binned entry rep- 
resents the results of a 30 minute collection cycle with nomi- 
nally identical high voltage value. The solid line represents a 
fitted Gaussian. 



D. Investigation of systematic errors 

We searched for a wide variety of potential systematic 
errors in our data in two basic ways. First, when two 
substantial subsets of the data could be identified, corre- 
sponding to two distinct experimental configurations, we 
compared the results of these subsets statistically. Sec- 
ond, in some situations, we could search for correlations 
in our Stark shift constant values as a function of some 
experimental parameter. In the first category, we com- 
pared data taken for up vs. down laser scans, for If 
vs. 2f demodulation, for F = 4 — > 5 vs. = 5 — > 4 
hyperfine transitions, for horizontal vs. vertical laser po- 
larization, and for HV on/off vs. off/on sequencing - 
none of which should influence our Stark shift constant. 
Fig. [To] summarizes these comparisons. Most of these 
searches yielded no statistically resolved systematic dif- 
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ferences. Note that the entire vertical extent of Fig 
represents 0.5% of our mean value. For the cases of laser 
sweep direction and laser polarization, where we see small 
{'^ 2 combined standard deviation) differences between 
subsets, we have included small systematic error contri- 
butions (see Table 1). 

We also searched for correlations in our data of the 
Stark shift constant with various experimental variables. 
Two examples of this are shown in Fig. [TT] At the top, 
we plot our measured values of kg versus electric field 






Full FM Fits 


V 


Polynomial Fits 


■ 


Up Scans 


□ 


Down Scans 


• 


Off ^ On 


o 


On ^ Off 


♦ 


1f demod. 





2f demod. 


▲ 


F=5 F'=4 


A 


F=4 ^ F'=5 


► 


Horz. Pol. 




Vert. Pol- 



FIG. 10. A summary of data subset comparisons to investi- 
gate potential systematic errors. The comparisons, in order, 
are: fit results based on complete FM line shape analysis vs. 
those derived from polynomial or Lorenztian fits to the central 
portion of the spectrum; direction of laser sweep; sequence of 
HV switching; RF lock-in demodulation frequency; hyperfine 
transition; and laser polarization. 
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FIG. 11. Correlation plots: fcg versus the square of the electric 
field (top), and ks versus laser power (below). 



(squared). We see no statistically resolved slope here, 
but we do observe a scatter in values slightly larger than 
expected based purely on statistical variation (reduced 
~ 2). For this reason we have included a small addi- 
tional systematic error contribution associated with the 
variation of kg with electric field in our final table of un- 
certainties. At the bottom of Fig. [Tl]we plot kg vs. laser 
power. One could worry that optical pumping or optical 
saturation effects or even AC Stark effects could affect 
our line shapes or line centers. While it is difficult to 
see how this could affect the HV on/off scans differently 
causing a systematic Stark shift error, it is nevertheless 
encouraging to see that over a factor of six variation in 
laser power we see no observable correlation. Finally, the 
list of error contributions in Table 1 includes a contri- 
bution from the overall uncertainty in the determination 
of the electric field (squared). This is dominated by our 
0.1% uncertainty in the plate separation. The other field- 
related uncertainties, from our high voltage divider and 
Keithley voltmeter, are an order of magnitude smaller. 
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We thus include a 0.2% contribution to the final error 
from the uncertainty in determination. As summa- 
rized in Table |Tj the quadrature sum of all errors, statis- 
tical and systematic, yields a 0.3% final uncertainty. 

TABLE I. Contributions to the final uncertainty. All entries 
are in units of kHz/(kV/cm)^. 





Sources of Error 


0.05 


Statistical 


0.05 


Electric Field Dependence 


0.14 


Laser scan direction 


0.20 


Laser polarization 


0.25 


E-Field Calibration 


0.33 


Quadrature Sum of Errors 



V. DISCUSSION OF RESULTS AND 
COMPARISON WITH THEORY 

Our final result for the Stark shift constant is: ks — 
-122.92(33) kHz/(kV/cm)2. Converting this result to 
scalar polarizability (in atomic units), we find that 
Aao = ao(65i/2) - ao(5Pi/2) = 1000.2 ± 2.7 ag. This 
result is in agreement with, but a factor of 30 more 
precise than, a measurement by Fowler and YellinJU] 
in which they found Aao — 944(73) Oq. As discussed 
in [IS], a new theory calculation of this quantity has 
just been completed. The result of this calculation gives 
Aao = 995(20) Cq, in excellent agreement with our ex- 
perimental result. 

The infinite sums in the theoretical calculations for 
the initial and final state polarizabilities are dominated 
(nearly 95% of the total contribution) by terms involv- 
ing mixing of the 6S1/2 state and the nearby 6Pi/2,3/2 
states. Because of this, we can combine the theoretical 
expression with our experimental result to infer highly 
accurate values for the matrix elements corresponding to 
these two transitions (and hence, the 6P excited state 
lifetimes). Using summations as in Eq. ([l]) as a start- 
ing point, we recast the expression for the polarizability 
difference to isolate the two dominant terms as follows: 

Aao^B-S + C, (4) 

where 5' = (6Pi/2||i:'||65)2 (the El line strength), and 
we define 



3\E{6Pi/2)- E{6S) E{6P3/2) - E{6S) J'^ ' 

in Eq. ^ is the ratio of the 65 - 6P3/2 to 65* - 
6P1/2 El line strength. In the present case, E? has the 
theoretical value [15] of 1.949(2), allowing us to compute 
a precise value for B. The constant C which accounts for 



all terms in the sums beyond the 65 — 6P mixing terms, 
has the theoretical value —60(8) ag in this case. Inserting 
our experimental value into the left hand side of Eq. Q , 
we can then infer a value for the line strength, S, with 
corresponding uncertainty (largely due to the theoretical 
error in C). Finally, using the known energy splittings, 
we can compute a value for the decay rate of the two 6P 
states, using the relation f?! 

_ 2.02613 X IQi^ Sab _i 

Aab- (2]:tt)' ' 

where the line strength. Sab, is expressed in atomic units, 
and the transition wavelength is in Angstroms. In this 
case, the 6P excited state lifetimes are given by Ta = 
1/Aab, since there is only one El decay channel from 
each of these states. We find 

r(6P3/2) = 58.17(45) ns 
r(6Pi/2) = 63.77(50) ns. 

These values, with uncertainties below 1%, are an order- 
of-magnitude more precise than earlier measurements de- 
rived from pulsed laser spectroscopy [H], and agree very 
well with recent theoretical estimates |4]. 



VI. CONCLUSIONS AND FUTURE WORK 

We have measured the scalar polarizability for the 
ground state 410 nm transition in atomic indium, improv- 
ing the experimental precision by more than an order of 
magnitude. This work has spurred a new set of atomic 
theory calculations whose results agrees very well with 
our measurement. We are currently introducing a second 
(infrared) laser system to our experimental setup. Using 
two-step excitation, we will then probe the 65i/2 — 6P3/2 
transition in our atomic beam with a diode laser centered 
at 1291 nm. We recently completed an indium vapor 
cell measurement using this exact optical arrangement fS] . 
We plan to use frequency modulation spectroscopy again 
to probe the very weak IR absorption signal. The Stark 
shift for this excited state transition is expected to be 
significantly larger, and will exhibit both scalar and ten- 
sor components, given the J = 3/2 nature of the excited 
state. Since the 6P3/2 state polarizability is strongly af- 
fected by the presence of nearby D-states, this presents 
a new challenge to atomic theory. Finally, having devel- 
oped this dual modulation scheme, effective for studying 
weak transitions in the atomic beam environment, we 
plan to study directly the thallium 6P1/2 — 6P3/2 Ml 
transition in the atomic beam. Here, because the va- 
por pressure of thallium is much larger, a similar atomic 
beam source can produce a much denser downstream 
beam, which would help to compensate for the intrin- 
sically small transition strength of the magnetic dipole 
line. 
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